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1  Objectives  and  results  of  the  research 

The  goal  of  the  project  “Modular  Paradigm  for  Scalable  Quantum  Information”  was  to  address  some  of 
the  challenges  facing  the  field  of  quantum  information  science  (QIS).  The  requirements  to  reliably  control 
a  scalable  quantum  system  while  staving  off  decoherence  pose  a  contradiction,  as  fast  control  implies  a 
strong  coupling  to  a  controlling  (external)  system,  but  this  entails  an  undesired  interaction  with  the  environ¬ 
ment,  leading  to  decoherence.  We  thus  studied  a  novel  paradigm  for  QIS,  given  by  smaller  building  blocks 
(composed  of  a  quantum  register  and  an  actuator)  that  are  connected  by  spin  wires. 

We  obtained  significant  results  in  both  these  tasks.  Over  the  course  of  the  project  we  developed  novel  control 
techniques,  such  as  Hamiltonian  engineering  for  many-body  systems  [1],  and  hybrid  control  combining 
classical  controllers  with  quantum  actuators  [2,  3].  A  significant  result  is  the  first  implementation  of  a 
feedback  control  algorithm  with  an  electronic  spin  qubit  in  diamond  [4].  We  experimentally  demonstrated 
protection  against  naturally  occurring  dephasing  noise,  extending  the  electronic  spin  qubit  coherence  time 
by  three  orders  of  magnitude.  This  is  a  critical  step  toward  building  fault-tolerant  quantum  devices  and 
it  will  also  contribute  to  the  wider  application  of  quantum  feedback.  In  addition,  we  explored  means  to 
connect  distributed  registers  via  spin  wires ,  by  performing  the  first  demonstration  of  quantum  information 
transport  with  a  mixed-state  spin  wire  [5].  We  further  explored  the  effects  of  noise  and  disorder  on  quantum 
information  transport,  elucidating  some  complex  dynamics  [6]  and  the  emergence  of  localization. 

2  Control  and  coherence  of  the  quantum  actuator/register  system 

The  first  activity  of  the  project  was  to  develop  control  strategies  for  the  actuator/register  system  and  to 
implement  them  experimentally. 

While  semiclassical  techniques  are  adequate  for  simple  manipulation  of  spins,  such  as  traditional  NMR 
studies  of  bulk  systems,  the  key  challenge  for  modern  quantum  engineering  is  to  achieve  near-perfect  con¬ 
trol  fidelity  in  the  presence  of  strong  coupling  to  the  environment.  To  meet  this  challenge,  we  moved  beyond 
the  control  paradigms  used  in  the  past,  devising  control  techniques  beyond  the  usual  rotating  wave  approxi¬ 
mation  [7]  and  combining  coherent  control  and  dissipation  for  state  preparation  [8]. 

An  innovative  theme  of  this  research  program  has  been  to  exploit  quantum  systems  themselves  as  controllers. 
Using  a  quantum  actuator  can  allow  faster  control;  it  avoids  additional  noise  sources  linked  to  a  classical 
apparatus;  and  it  can  achieve  a  higher  degree  of  spatial  resolution. 

We  thus  studied  time-optimal  control  via  a  quantum  actuator  [2]  and  investigated  when  this  strategy  can 
beat  driving  by  a  classical  field  [9].  Significantly,  we  demonstrated  experimentally  that  combining  quantum 
and  classical  controllers  can  achieve  even  faster  control  [3].  A  quantum  controller  is  often  also  the  best 
option  to  implement  feedback  control,  a  particularly  efficient  strategy  to  protect  qubits  from  detrimental 
noise  while  implementing  desired  operations.  Open-loop  control  such  as  dynamical  decoupling  suppresses 
slowly  varying  phase  noise  and  enhances  the  qubit  coherence  time.  However,  it  is  ineffective  against  fast- 
varying  noise  and  often  lacks  compatibility  with  quantum  gates.  Feedback-based  control  could  overcome 
these  difficulties  and  is  indeed  the  method  of  choice  for  the  control  of  classical  systems.  Nonetheless,  as 
feedback  usually  entails  destructive  measurements  that  complicate  the  qubit  dynamics,  until  our  work  the 
implementation  of  feedback  protocols  had  been  limited  mainly  to  quantum  optics.  Using  the  NV  nitrogen 
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nuclear  spin  as  an  ancillary  qubit,  we  experimentally  demonstrated  protection  against  naturally  occurring 
dephasing  noise,  extending  the  electronic  spin  qubit  coherence  time  by  three  orders  of  magnitude  [4].  This 
feedback  scheme  is  compatible  with  applications  of  gates  on  the  protected  qubit,  as  we  demonstrated  by 
performing  a  protected  NOT  gate.  I  expect  that  this  demonstration  of  the  power  of  feedback  control  for 
the  precise  manipulation  of  quantum  systems  will  contribute  to  its  widespread  application,  analogous  to  its 
success  in  classical  control. 

Among  our  results, 

•  We  found  the  shape  of  the  time-optimal  control  strategy  to  engineer  desired  qubit  gates  via  the  actuator 
control.  Using  algebraic  methods  we  were  able  to  find  the  time-optimal  control  solutions  for  the  most 
general  SU(2)  synthesis  problem. 

•  We  derived  new  bounds  on  the  time  and  number  of  switches  required  to  engineer  desired  qubit  gates  via 
the  actuator  control  that  are  more  general  and  stricter  than  previously  known  bounds. 

•  We  developed  an  efficient  numerical  method  to  find  the  time-optimal  control  sequence  and  applied  it  to 
find  gates  for  the  coupled  NV-nuclear  spin  system. 

•  We  devised  a  method  to  implement  time-optimal  bang-bang  control  in  the  practical  situations  where  the 
experimental  control  doesn’t  have  infinite  bandwidth  and  quantified  the  ensuing  loss  of  fidelity. 

•  We  characterized  experimentally  the  nitrogen  nuclear  spin  Hamiltonian,  devising  a  method  to  measure 
the  full  hyperfine  tensor.  This  measurement  had  not  been  possible  until  now  for  single  NV  centers. 

•  We  exploited  the  control  achieved  on  the  nitrogen  nuclear  spin  to  implement  a  feedback  algorithm  to 
preserve  the  coherence  of  the  NV  electronic  spin 

Details  of  these  activities  are  described  below. 

2.1  Optimal  control  of  a  qubit  by  a  quantum  actuator 

A  central  goal  of  the  project  was  to  devise  and  demonstrate  quantum  control  techniques  where  the  control¬ 
ling  apparatus  itself  is  a  quantum  system.  Using  a  quantum  controller  can  allow  faster  dynamics;  it  avoids 
additional  noise  sources  linked  to  a  classical  apparatus;  and  it  can  achieve  a  higher  degree  of  spatial  reso¬ 
lution.  We  focused  our  attention  to  the  control  of  well-isolated  single  spin  qubits  by  a  quantum  controller 
(itself  a  spin  qubit),  which  is  strongly  coupled  to  external  driving  fields.  The  goal  is  the  generation  of  desired 
unitaries  (or  gates)  to  control  the  system  for  any  initial  state.  The  Hamiltonian  of  the  system  can  be  written 
as 

Uq  =  uLaz  +  SazA  ■  a,  (1) 

where  S %  refers  to  the  actuator  spin  and  a  to  the  nuclear  spin  in  the  register.  A  simple  control  strategy  is 
to  let  the  actuator  alternate  between  two  possible  states  by  applying  7r-pulses  at  designated  times  t Then 
the  Hamiltonian  (1)  alternates  between  two  non-commuting  operators  z  =  and  v  =  ujl<jz  +  A  •  a. 
Switching  between  the  actuator  eigenstates  is  enough  to  achieve  full  controllability  of  the  qubit,  as  long  as 
the  coupling  is  anisotropic  [10,  11].  In  the  absence  of  direct  driving  of  the  qubit,  Pontryagin’s  minimum 
principle  proves  that  this  bang-bang  control  achieves  time-optimality  [12-14]  and  one  needs  to  find  a  series 
of  switching  times  {£&}  that  give  the  desired  evolution. 

More  generally,  we  aim  at  solving  the  problem  of  generating  a  transformation  in  SU(2)  [or  SO(3)]  when  the 
available  control  is  given  by  rotations  around  two  non-parallel  axes,  x  and  v  =  ft  (cos  ax  +  sin  ay),  with 
ft  >  1,  with  clockwise  ( t &  >  0)  rotations  only  or  both  clockwise  and  anti-clockwise  rotations  ( t &  ^  0). 
Thus,  any  allowed  control  sequence  can  be  written  as 

U  =  x(h)  V(t2)  X(tx>k)  V(tv,k+ 1)  . . . ,  with  X(t)  =  V(t)  =  e~i3^/2  (2) 

We  used  an  algebraic  framework  to  study  this  time-optimal  synthesis  problem,  which  is  common  in  both 
quantum  and  classical  mechanics.  The  method  we  used  bypasses  usual  control-theoretical  techniques,  and 
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easily  imposes  necessary  conditions  on  time-optimal  sequences.  We  thus  derive  the  fact  that  time-optimal 
sequences  are  solely  parametrized  by  three  rotation  angles  and  a  finite  combination  of  signs  [2]. 


Case  t  >  0.  Time-optimal  sequences  only  depend  on  four  parameters,  namely  the  total  number  of  rotations 
n  <  oo,  the  outer  angles  t\An  and  the  internal  angle  tx  (or  tv ).  In  this  case,  the  internal  angles  are  related 
by 


f  tx\  hi  —  cos  (a) 
\2  J  1  —  hicos(a) 


(14) 


Case  t  ^  0.  It  holds  that  either  n  <  5,  or  n  oo.  Time-optimal  sequences  only  depend  on  four  parameters, 
as  above,  and  a  few  combinations  of  signs  (only  four  signs  need  to  be  specified).  Internal  angles  are  related 
by 

tan(|)=±tlm(|)i  (18) 

in  finite  sequences;  and,  in  infinite  sequences 


tan 


tan 


n  —  cos(o') 

1  —  ttcos(a') 

i  hi  +  cos  (a) 

'  1  +  hicos(a) 


(22) 

(23) 


Furthermore,  we  were  able  to  derive  general  bounds  on  the  times  as  a  function  of  the  angle  between  the  axes 
and  the  relative  rotation  speed  of  each  control.  Results  are  substantially  different  whether  both  clockwise 
and  counterclockwise  rotations  about  the  given  axes  are  allowed,  or  only  clockwise  rotations.  In  the  first 
case,  we  proved  that  any  finite  time-optimal  sequence  is  composed  at  most  of  five  control  concatenations, 
while  for  the  more  restrictive  case,  we  found  scaling  laws  on  the  maximum  length  of  any  finite  time-optimal 
sequence;  bounds  for  both  cases  are  stricter  than  previously  published  ones.  These  bounds  severely  constrain 
the  structure  of  time-optimal  sequences,  allowing  for  a  simple  numerical  search  of  the  time-optimal  solution. 
We  note  that  besides  the  application  to  our  present  quantum  control,  the  theory  is  much  more  general  and 
can  be  applied  even  to  classical  rotations,  for  example  for  robotics  applications. 

Our  theoretical  results  can  be  immediately  used  to  the  problem  of  controlling  nuclear  spins  coupled  to  the 
NV  center  in  diamond.  As  a  concrete  example,  we  consider  a  single  NV  center  electronic  spin  S  =  1 
coupled  to  a  13C  nuclear  spin  7=1/2.  Their  Hamiltonian  is 


H  =  A S2Z  +  7 eB0Sz  +  1CB0IZ  +  S  •  A  •  I , 


(3) 


where  A  =  2.87GHz  is  the  NV  zero-field  splitting;  7e  «  2.8MHz/G,  7 c  ~  lkHz/G  are,  respectively,  the 
gyromagnetic  ratios  of  the  electron  and  nuclear  spins;  7>o  is  a  static  magnetic  field  along  the  NV  5-axis;  and 
A  is  the  hyperfine  tensor.  The  NV  spin  triplet  can  be  reduced  to  an  effective  two-level  system  by  driving  the 
system  on  resonance  with  a  transition  between  two  Sz  eigenstates  (e.g.  | ms  =  0)  \ms  =  +1)),  while  the 
third  eigenstate  (e.g.  \ms  =  —1))  can  be  neglected.  Then  the  Hamiltonian  can  be  rewritten  in  the  electronic 
spin  rotating  frame  as 

U  =  u)QIz  +  SZAZ  •  1=  \0)(0\ujoIz  +  I±l)(±l|  (lj0Iz  ±  Az  •  I)  ,  (4) 


where  ujq  =  jcBq  (that  we  assume  >  0).  The  contact  and  dipolar  contributions  [15]  to  the  hyperfine 
coupling  A  can  be  described  by  a  longitudinal  component  Ay  and  a  transverse  component  A±,  that  we  will 
take  without  loss  of  generality  along  the  x  direction.  The  nuclear  spin  thus  rotates  around  two  distinct  axes, 
depending  on  the  electronic  spin  manifold.  Then,  a  simple  strategy  for  the  indirect  control  of  the  nuclear 
spin  is  to  induce  alternating  rotations  by  flipping  the  electronic  spin  state  with  (fast)  7r-pulses.  We  define  the 
axes  and  rotation  speeds  in  the  two  manifolds  as 

W0  =  1C B0  =  KW±1,  W±1  =  y/(w0±yl||)2  +  ^j_  ,  ^ 

vo  =  5  ,  v\  =  5cos(a)  +  xsin(a)  , 
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Fig.  1:  Comparison  of  gate  time,  Left:  Case  k  <  cos(a),  occurring  for  a  13C  at  a  distance  of^  2.92A  from  the  NV 
center,  with  an  external  magnetic  field  Bq  «  500G  aligned  with  the  z  axis.  We  plot  the  simulated  actuator  implemen¬ 
tation  time  (blue  circles-left  axis)  of  the  unitaries  X($)  (left)  and  Y($)  (right)  and  the  corresponding  sequence  lengths 
(red  crosses-right  axis).  For  comparison,  we  plot  the  time  required  with  direct  driving  (green  lines)  with  bare  Rabi 
frequencies  20  and  100kHz,  when  the  electronic  spin  in  state  |  —  1)  (left),  thus  maximizing  the  enhancement  factor,  or 
|0)  (right).  Note  that  the  direct-driving  time  for  $  >  it  depends  on  whether  the  driving  phase  can  be  inverted  (dashed 
line)  or  not  (solid  line).  Comparison  of  gate  time,  Right:  Case  k  >  cos(o),  occurring  for  a  13C  at  a  distance  of 
«  4.3lA  from  the  NV  center.  Note  that  virtual  transition  of  the  electronic  spin  in  the  ms  =  0  manifold  result  in  a 
decrease  of  the  effective  Rabi  frequency,  thus  making  direct  driving  in  that  manifold  unfavorable. 


Wlth  tan(a)  =  — ~rj~  >  K  =  ~  •  (6) 

(jJ 0  ±  A||  UJ±\ 

If  the  NV  electronic  spin  is  initially  in  the  |0)  state,  applying  7r-pulses  at  times  T k  gives  the  nuclear  spin 
evolution: 

JJ  —  e-i<PnVi-v  m  '  '  e-i<PkV0-° e-iVk- 1^1-°  m  m  m  e-i<PiV0-v^ 

where  j>°k^  =  (T&  —  for  odd  (even)  fc,  and  a  are  the  Pauli  matrices. 

An  alternative  strategy  for  qubit  control  is  to  use  classical  driving  fields.  Resonant  driving  along  a  desired 
rotation  axis  achieves  time-optimal  steering  of  the  qubit  in  the  xy  plane  [12,  16].  Even  when  the  direct 
driving  of  the  qubit  is  slow,  the  rate  might  be  increased  by  virtual  transition  of  the  actuator.  This  is  the  case 
for  nuclear  spins:  while  their  coupling  to  an  external  driving  field  is  weak,  indirect  forbidden  transitions 
mediated  by  the  electronic  spin  can  considerably  enhance  the  driving  strength  [17-19].  This  nuclear  Rabi 
enhancement  depends  on  the  state  of  the  electronic  spin.  The  effective  Rabi  frequency  Q  for  an  isolated 
nuclear  spin,  hence,  is  modified  from  its  bare  value  Q  by  the  enhancement  factors  Co,±i>  corresponding  to 
the  electronic  spin  states  |0),  |  ±1).  The  enhancement  is  proportional  to  the  ratio  of  the  qubit  and  actuator 
coupling  to  the  external  field.  For  nuclear  and  electronic  spins  considered  here,  je/ln  ~  2600  and  the 
effective  Rabi  frequencies  fA  =  ( 1  +  Q)Q  can  be  much  larger  than  the  bare  frequency. 

We  assume  «  100kHz  as  an  upper-limit  on  realistic  bare  nuclear  Rabi  frequencies  by  considering  data 
in  [20],  where  the  13  C  considered  was  only  weakly  coupled  and  thus  no  Rabi  enhancement  was  present. 
To  achieve  this  strong  driving,  a  dedicated  microfabricated  coil  was  necessary  [21].  Rabi  frequencies  « 
20kHz  are,  in  our  experience,  in  the  upper  achievable  range  with  modest  amplifiers  and  a  simple  wire  to 
deliver  the  rf  field. 

Both  regimes  of  k  ^  cos (cr)  for  the  time-optimal  solutions  can  be  explored  in  the  NV  center  system  by 
considering  the  coupling  to  13  C  at  different  distances  from  the  NV  defect  [22-24].  The  hyperfine  tensors 
for  13  C  located  up  to^  8 A  away  from  the  NV  center  were  estimated  using  density  functional  theory  [25]. 
In  what  follows,  we  numerically  compare  the  performance  of  the  proposed  control  method  against  direct 
driving  under  diverse  experimental  conditions  and  for  a  number  of  distinct  nuclear  spins.  For  the  NV  center 
system  the  dependence  on  the  hyperfine  parameters  of  both  the  actuator  scheme  time  and  the  direct  driving 
strength  yields  a  broad  variation  of  results  for  both  close-by  and  more  far  away  nuclei;  while  a  trend  toward 
longer  times  for  the  actuator  scheme  vs.  direct  driving  is  apparent  as  the  distance  from  the  NV  center 
increases,  the  large  variations  indicate  that  the  best  scheme  should  be  evaluated  for  individual  nuclear  spins. 
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Fig.  2:  Left:  energy  levels  of  the  reduced  NV-14N  spin  system,  showing  the  transitions  that  are  mixed  by  the  transverse 
hyperfine  coupling.  Center:  14N  Rabi  Frequency  in  the  three  NV  manifold  (Red,  solid  line  ms  =  0.  Black,  dashed 
line,  ms  =  —1.  Gray,  dotted  line  ms  =  +1)  as  a  function  of  the  magnetic.  The  filled  symbols  correspond  to 
the  experimental  data,  which  matches  closely  the  theoretical  prediction.  The  effective  Rabi  frequencies  increase 
rapidly  with  the  field,  exceeding  1MHz  when  close  to  ground  state  level  anti-crossing.  The  enhancement  allows 
fast  manipulation  of  the  nuclear  spin  even  when  the  bare  Rabi  field  is  only  B\  «  3G.  The  theoretical  prediction  is 
confirmed  by  simulations  (open  symbols)  of  the  spin  dynamics.  Right:  14N  Rabi  oscillations  at  B  =  450G  in  the  three 
NV  manifold  (Red,  solid  line  ms  =  0.  Black,  dashed  line,  ms  =  —1.  Gray,  dotted  line  ms  =  +1).  Here  the  dots  are 
the  experimental  results,  while  the  lines  are  fits  to  cosine  oscillations.  The  different  baseline  of  the  ms  =  —  1  curve  is 
due  to  small  differences  in  the  fluorescence  emission  of  different  nuclear  manifolds. 


Indirect  control  of  qubits  by  a  quantum  actuator  is  an  attractive  strategy  in  many  situations  when  the  qubits 
couple  weakly  to  external  fields,  but  interact  more  strongly  to  another  quantum  system. 

An  interesting  extension  of  our  results  would  be  to  simultaneously  control  two  or  more  qubits  by  the  same 
quantum  actuator.  While  this  is  possible,  provided  the  qubits  are  coupled  with  different  strengths  [10], 
it  becomes  more  difficult  to  find  time-optimal  solutions  except  for  particular  tasks  (such  as  state-to- state 
transformations  [26])  or  geometries  [27,  28].  Still,  even  when  the  goal  is  to  control  a  larger  number  of 
qubits,  our  results  can  guide  the  experimentalist’s  choice  between  direct  driving  and  the  actuator  control,  for 
which  these  results  give  an  upper  bound. 

2.2  Hybrid  control:  direct  driving  enhanced  by  a  quantum  actuator 

When  both  direct  driving  and  a  quantum  controller  are  available,  one  can  take  advantage  of  both  to  achieve 
faster  driving  of  the  qubit.  In  particular,  we  found  that  the  quantum  actuator  can  mediate  virtual  transitions 
that  enhance  the  direct  driving  strength.  In  turn,  since  the  enhancement  depends  on  the  actuator/register  cou¬ 
pling  strength,  the  enhancement  can  be  used  to  precisely  measure  the  coupling  itself  [3].  Precise  knowledge 
of  a  system’s  Hamiltonian  is  a  pre-requisite  for  its  accurate  control. 

The  NV  ground  state  is  a  two-spin  system  given  by  the  electronic  spin  of  the  NV  center  ( S  =  1)  and 
the  nuclear  spin  (/  =  1)  of  the  substitutional  14N  adjacent  to  the  vacancy  that  comprise  the  defect.  In 
the  experiments,  we  are  only  interested  in  two  of  the  nuclear  spin  levels  (m/  =  +1,0)  that  we  drive  on- 
resonance,  while  the  third  level  can  be  neglected.  Then,  the  Hamiltonian  of  the  reduced  system  is  given  by 
7+  =  7+||  +  T~L±,  where  the  secular,  Hy,  and  nonsecular,  1~L±,  terms  are: 

H\\  =  ^S2z  +  {leBz  +  ^-)Sz  +  {Q  +  lnBz)Iz  +  A\\SzIz, 

n±  =  V2A±(SXIX  +  Syly).  (8) 

Here  S  and  /  are  the  electron  spin-1  and  nuclear  spin- 1/2  operator  respectively,  A  =  2.87  GHz  is  the  zero- 
field  splitting  and  Q  =  —4.945  MHz  [29]  the  nuclear  quadrupolar  interaction.  The  NV  spin  is  coupled  to 
the  nuclear  spin  by  a  hyperfine  interaction  with  a  longitudinal  component  Ay  =  —2.162  MHz  [29]  and  a 
transverse  component  A±  which  we  want  to  estimate.  A  magnetic  field  Bz  is  applied  along  the  NV  crystal 
axis  [1 1 1]  to  lift  the  degeneracy  of  the  ms  =  ±1  level,  yielding  the  electron  and  nuclear  Zeeman  frequencies 
7 eBz  and  ynBz  where  ye  =  2.8  MHz/G  and  yn  =  0.308  kHz/G. 

Let  |ras,raj)  be  eigenstates  of  7/y.  The  transverse  coupling  A±  mixes  states  connected  via  zero-quantum 
(ZQ)  transitions,  |+1,0)  «+  |0, 1)  and  |0,0)  o-  |  — 1,1).  Thus,  diagonalization  of  the  total  Hamiltonian 


DISTRIBUTION  A:  Distributicgi  approved  for  public  release. 


can  be  achieved  by  rotating  the  two  ZQ  subspaces  with  a  unitary  transformation  Uzq  =  e~l^y  0  +cry°+\ 
where  the  rotation  operators  are  given  by  =  z(|+l,0)(0, 1|  —  |0, 1)(+1,0|);  a~  =  z(|0, 0)(— 1, 1|  — 
|  — 1, 1)  (0,  0|)  and  the  rotation  angles  are 


tan(20+) 

tan(20-) 


2A± 

A  +  7 eBz  -  7 nBz  -  Q 
-2A± 

A  —  jeBz  —  An  +  7 nBz  +  Q 


(9) 

(10) 


Because  of  this  level  mixing,  a  field  on  resonance  with  the  nuclear  spin  transition  also  drives  electronic 
transitions.  Although  the  electronic  spin  state  is  unchanged  to  first  order,  as  long  as  the  mixing  is  small,  the 
forbidden  transitions  result  in  an  enhancement  of  the  nuclear  state  driving  frequency,  as  we  explain  below. 
When  applying  a  radio  frequency  (RF)  field  to  drive  the  nuclear  spin,  the  interaction  Hamiltonian  of  the 
NV-14N  system  with  the  RF  field  is: 


=  2Bi  cos(ujt)(-feSx  +  V2ZnIx), 


(11) 


where  B\  is  the  RF  field  strength.  The  Hamiltonian  can  be  simplified  by  going  into  a  rotating  picture  at  the 
RF  frequency  uo  and  applying  the  rotating  wave  approximation  (RWA),  to  obtain  Br f  =  B\  {^eSx+V2^nIx). 
We  note  that  since  we  might  have  7e£>i  cc,  effects  from  the  counter-rotating  fields,  such  as  Bloch-Siegert 
shifts  of  the  electronic  energies,  might  be  present.  These  effects  were  however  negligible  at  the  fields  and 
Rabi  strengths  used  in  the  experiments.  Transforming  Br f  with  the  unitary  U  and  denoting  the  states  and 
operators  in  the  new  frame  by  a  hat,  we  obtain  BY{  =  UzqBy{ (£)^zq  =  Bn  +  Be,  with 


'Hn  =  'fnBi  (ai|i)(i|e+a0|6)(6|e+a_i|-i)(-i|e)  Ix 

Here  ams  denote  the  enhancement  factors  in  each  manifold  of  the  NV  spin: 


cr+i  ~  1  + 

CX  o  ~  1  — 


+ 

( 1  I  ^  1  -)- 


7e 

A± 

In 

A  +  r)eBz  rynBz 

Q’ 

7e 

(  A ± 

In 

V  A  +  7 eBz  -  7 nBz  - 

-Q 
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+  Q )’ 

7e 

A± 

Tn  ^  1 e^z  Ah  +  ZnBz  +  Q 


(12) 


(13) 


(14) 

(15) 


where  we  show  expressions  exact  up  to  the  first  order  in  6±  (see  [3]  for  the  exact  expressions).  The  Hamil¬ 
tonian  Be  can  be  neglected  since  electronic  spin  transitions  are  far  off-resonance. 

To  confirm  the  expected  dependence  of  the  Rabi  enhancement  factors  on  the  external  magnetic  field  and  the 
NV  state,  we  measured  the  Rabi  oscillations  at  the  three  electronic  spin  manifolds  with  varying  magnetic 
field  Bz.  As  shown  in  Fig.  (2),  the  measured  Rabi  frequencies  match  well  with  the  theoretical  model.  It 
is  worth  noting  that  contrary  to  the  static  pseudo-nuclear  Zeeman  effect  [30],  there  is  a  large  enhancement 
(rro  ~  16,  a±i  «  —9)  even  at  zero  field.  Also,  close  to  the  ground  state  avoided  crossing  (B  «  0.1  T) 
the  enhancement  can  become  very  large,  exceeding  100.  The  validity  of  our  approximation  in  this  regime 
can  be  confirmed  by  numerical  simulations,  based  on  Trotter  expansion  and  Floquet  theory.  Thanks  to  the 
strong  dependence  of  the  enhancement  factors  from  the  transverse  hyperfine  coupling,  we  can  determine 
A_\_  from  measurement  of  the  14N  Rabi  oscillations  with  higher  precision  than  ever  done  previously,  even  in 
ensemble  measurements. 
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2.3  Quantum  Feedback 


Engineering  desired  operations  on  qubits  subjected  to  the  deleterious  effects  of  their  environment  is  a  critical 
task  in  quantum  information  processing,  quantum  simulation  and  sensing.  The  most  common  approach  is 
to  rely  on  open-loop  quantum  control  techniques,  including  optimal  control  algorithms,  Lyapunov  design 
and  Hamiltonian  engineering  [1].  An  alternative  strategy,  inspired  by  the  success  of  classical  control,  is 
feedback  control.  Because  of  the  complications  introduced  by  quantum  measurement,  closed-loop  control 
is  less  pervasive  in  the  quantum  settings  and  its  experimental  implementations  have  been  mainly  limited  to 
quantum  optics  experiments.  We  implemented  a  feedback  control  algorithm  with  a  solid-state  spin  qubit 
system  associated  with  the  Nitrogen  Vacancy  (NV)  centre  in  diamond,  using  coherent  feedback  to  over¬ 
come  limitations  of  measurement-based  feedback  (see  Fig.  7),  and  show  that  it  can  protect  the  qubit  against 
intrinsic  dephasing  noise  for  milliseconds  (Fig.  4). 

In  coherent  feedback,  the  quantum  system  is  connected  to  an  auxiliary  quantum  controller  (ancilla)  that  ac¬ 
quires  information  about  the  system’s  output  state  (by  an  entangling  operation)  and  performs  an  appropriate 
feedback  action  (by  a  conditional  gate).  In  contrast  to  open-loop  dynamical  decoupling  (DD)  techniques, 
feedback  control  can  protect  the  qubit  even  against  Markovian  noise  and  for  an  arbitrary  period  of  time  (lim¬ 
ited  only  by  the  ancilla  coherence  time),  while  allowing  gate  operations.  It  is  thus  more  closely  related  to 
Quantum  Error  Correction  schemes,  which  however  require  larger  and  increasing  qubit  overheads.  Increas¬ 
ing  the  number  of  fresh  ancillas  allows  protection  even  beyond  their  coherence  time  as  well  as  protection 
against  multi-axis  noise. 

We  further  evaluated  the  robustness  of  the  feedback  protocol,  which  could  be  applied  to  quantum  computa¬ 
tion  and  sensing,  by  exploring  an  interesting  tradeoff  between  information  gain  and  decoherence  protection, 
as  measurement  of  the  ancilla-qubit  correlation  after  the  feedback  algorithm  voids  the  protection,  even  if  the 
rest  of  the  dynamics  is  unchanged. 
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Fig. 3:  Feedback  algorithm  and  experimen¬ 
tal  implementation,  (a)  NV  centre  in  diamond 
and  (b)  relevant  energy  levels  of  the  spin  sys¬ 
tem,  showing  polarisation  processes  under  op¬ 
tical  illumination  [31].  (c)  Quantum  circuit. 

Hadamard  gates  prepare  and  read  out  a  super¬ 
position  state  of  the  qubit,  \(p)q  =  -F(|0)  + 

|1)).  Amid  entangling  gates  between  qubit  and 
ancilla,  the  qubit  is  subjected  to  noise  (and  pos¬ 
sibly  unitary  gates  U).  We  assume  the  ancilla 

is  not  affected  by  the  bath,  yielding  ta0Jqb(r),  with  Uqb(r)  the  qubit-bath  joint  evolution.  Given  a  dephasing  bath,  we 
set  the  entangling  gate  to  Uc  =  crx  (conditional-NOT  gate).  More  generally,  upon  undoing  the  entangling  operation, 
the  system  is  left  in  the  state  |^(r))  =  ^=(|0)a  K+  \ ipq,£h)  +  |l)a  K~  \ ipq,£b)),  with  K*1  =  Uqb  ±  UcUqhU\.  The 
entangling  gate  Uc  is  designed  such  that  K+  =  lq  ®  \b  and  ®  Xb  ,  where  xt  act  on  the  bath  only,  and  Uq 

on  the  qubit.  After  measuring  the  ancilla,  we  could  use  a  feedback  operation  Uq  to  restore  the  correct  qubit  state.  The 
ancilla  measurement  is  replaced  by  coherent  feedback  (shaded  region)  obtained  by  a  controlled-correction  gate  (here 
Uq  =  ct z  for  dephasing  noise).  The  final  state  of  the  combined  system  is  then  -^(\$)aXb  +  \1) a  X.b  )  \v) q  \€) b>  which 
reveals  how  the  qubit  is  now  decoupled  from  the  bath. 

(d)  Experimental  implementation.  The  laser  excitation  polarises  both  spins.  Black  sinusoidal  lines  refer  to  selective 
MW  pulses  acting  only  in  the  raj  =  1  manifold  (thus  mimicking  controlled-NOT  gates),  while  solid  bars  indicate 
non-selective  pulses.  The  RF  excitation  describes  selective  pulses  in  the  ms  =  0  manifold.  We  use  7t/2  rotations 
about  x  to  approximate  Hadamard  gates.  To  implement  a  nonselective  RF  7t/2  gate  on  the  nuclear  spin  we  embed  a 
nonselective  MW  i r  pulse  into  two  consecutive  RF  7t/2  pulses.  The  controlled-correction  gate  is  implemented  by  free 
evolution  ( tz )  under  hyperfine  coupling. 
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Fig.  4:  Experimental  demonstration  of  the  feedback-based  protection  algorithm,  (a)  The  signal  -normalised 
Photoluminescence  (PL)  intensity-  oscillates  at  the  hyperfine  coupling  frequency,  A  =  —2.15  MHz.  The  initial 
coherent  superposition  state  of  the  qubit  is  preserved  for  a  time  r  >  1  ms  at  390  G  (red  dots),  while  we  observe 
a  sharp  decrease  in  the  signal  amplitude  at  514  G  (black  squares),  where  correlations  between  the  qubit  and  ancilla 
states  are  partly  measured.  This  is  evident  in  the  lower  panels,  where  we  compare  the  fidelity  at  short  (b),  and  long 
(c)  times,  for  B  =390  G  and  514  G.  To  highlight  the  differences  while  taking  into  account  different  PL  intensities 
at  the  two  fields,  we  normalised  all  the  data  so  that  at  short  times  the  signal  has  the  same  (maximum)  contrast,  (d) 
Protected  NOT  gate  (red  circles).  Here  we  show  that  the  coherence  of  the  qubit  is  protected  for  a  time  longer  than  the 
dephasing  time,  r  >  T£ e,  even  when  a  NOT  gate  is  applied.  We  compare  the  dynamics  to  the  NOOP  dynamics  (black 
squares),  clearly  showing  that  the  NOT  gate  inverts  the  state  of  the  qubit,  as  signalled  by  the  out-of-phase  oscillations, 
(e)  Weak  measurement  of  the  ancilla:  normalised  PL  signal  after  a  protection  time  r  =  8/iS,  as  a  function  of  the  ancilla 
measurement  strength.  In  the  experiment,  we  vary  the  ancilla  measurement  strength  by  changing  the  angle  of  the  last 
controlled  phase  rotation  gate. 


The  success  of  the  feedback-based  protection  algorithm  rests  on  the  fact  that  the  increased  qubit  entropy, 
due  to  the  coupling  to  the  environment,  is  dumped  on  the  ancilla.  This  could  be  revealed  by  measuring 
the  state  of  the  ancilla,  which  would  yield  information  about  the  noise,  while  preserving  the  qubit  state.  If 
instead  information  about  the  ancilla  is  collected  by  a  correlated  qubit-ancilla  measurement,  the  protection 
fails  [32],  as  it  would  happen  in  a  measurement-based  feedback  for  an  imperfect  ancilla  readout.  We  thus 
investigated  this  tradeoff  between  the  protection  power  and  information  gain  on  the  ancilla.  We  transferred 
a  part  of  the  ancilla  entropy  back  to  the  qubit  by  employing  a  conditional  gate  that  maps  the  state  of  the 
ancilla  onto  the  qubit,  correlating  the  two  qubit  states.  By  changing  the  angle  of  rotation  of  the  conditional 
gate  we  can  vary  the  strength  of  the  ancilla  measurement,  from  a  weak  measurement  to  a  strong  one.  As 
more  information  about  the  ancilla  is  acquired,  the  fidelity  of  the  protection  gate  decreases  (Fig.  4b).  Note 
that  we  can  combine  the  conditional  gate  performing  the  weak  measurement  of  the  ancilla  state  with  the  last 
conditional  Z-gate  of  the  algorithm,  thus  in  practice  performing  a  conditional  phase-shift  gate.  The  qubit 
fidelity  is  maximised  when  the  conditional  gate  performs  the  required  7r-rotation  on  the  qubit,  whereas  it 
decreases  when  a  different  phase  rotation  is  employed. 

Our  results  provide  important  information  about  the  robustness  of  quantum  feedback. 
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3  Spin  wires  for  quantum  information  transfer 


The  development  of  small  quantum  registers  must  be  complemented  by  a  strategy  to  couple  them.  Thus, 
the  second  activity  of  the  program  explored  quantum  information  transfer  using  quantum  spin  wires,  which 
transport  information  via  the  free  evolution  under  their  mutual  interaction. 

In  this  project  we  theoretically  studied  conditions  for  perfect  quantum  state  transfer  and  we  investigated  an 
experimental  implementation  with  nuclear  magnetic  resonance  techniques  applied  to  apatite  crystals.  We 
further  studied  the  effects  of  decoherence  on  spin  transport  and  we  the  performance  of  decoupling  techniques 
in  counteracting  the  effects  of  a  spin  bath.  In  addition,  we  looked  at  different  experimental  systems  that  could 
provide  a  platform  for  quantum  information  transport. 

Achieved  results  include: 

•  We  performed  the  first  experimental  implementation  of  quantum  information  transfer  in  mixed- state  spin 
chains. 

•  We  studied  the  coherence  properties  of  highly  correlated  spin  states  under  the  effects  of  a  correlated 
noise.  We  found  that  the  restriction  to  one-dimensional  geometries  brings  both  a  complex,  non-Markovian 
dynamics  and  also  longer  coherence  times  that  in  3D  systems,  thus  pointing  to  advantages  to  be  found  in 
particular  geometries  for  larger  quantum  information  architectures. 

•  We  investigated  the  effect  of  dephasing  and  depolarizing  noise  on  quantum  state  transfer.  For  this  study, 
we  devised  multi-pulse  control  sequences  that  re-introduce  in  a  controlled  way  the  effects  of  a  spin  bath 
onto  the  spin  chains.  This  allows  us  to  engineer  different  spin-bath  effects 

•  We  studied  the  polarization  transfer  from  NV  centers  in  diamond  (that  can  be  polarized  optically)  to 
dark  spins  associated  with  single  Nitrogen  defects  (PI  centers)  in  diamond.  We  are  now  studying  the 
polarization  dynamics  in  the  PI  spin  system,  with  the  goal  of  assessing  its  potential  as  a  bus  for  quantum 
information  transfer. 

Details  of  these  activities  are  described  below. 

3.1  Experimental  implementation  of  quantum  information  transfer 

We  used  nuclear  spin  systems  in  apatite  crystals  as  a  test-bed  to  probe  quasi-one-dimensional  (ID)  dy¬ 
namics,  including  transport  and  decoherence,  and  in  particular  to  study  protocols  for  quantum  information 
transport.  The  crystal  structure  of  fluorapatite  [Cas^O^F,  FAp]  and  hydroxylapatite  [Ca5(P04)3(0H), 
HAp]  presents  a  favorable  geometry  where  19 F  or  1H  nuclear  spins  are  aligned  in  linear  chains  along  the 
crystal  c-axis  with  inter-spin  spacings  much  shorter  than  the  distance  to  other  parallel  chains.  The  system 
thus  acts  as  an  ensemble  of  linear  spin  chains  with  well  characterized  internal  Hamiltonian.  NMR  techniques 
enable  the  exploration  of  quantum  transport  even  in  the  absence  of  single-spin  addressing  and  readout.  We 
used  a  superconducting  magnet  (Oxford  instrument,  7  T,  wide-bore)  and  a  spectrometer  (Bruker  Avance300 
DMX)  to  control  the  nuclear  spin  system.  The  spectrometer  is  optimized  for  solid-state  systems,  achieving 
high-power  pulses  (300  W  amplifiers)  to  address  the  broad  linewidth  of  solid-state  crystals  as  well  as  fast 
acquisition. 

While  an  Hamiltonian  for  quantum  state  transfer  can  be  relatively  easily  obtained  with  multiple  quantum 
sequences,  it  is  more  challenging  to  initialize  and  readout  the  system  in  the  absence  of  single-spin  address¬ 
ability.  Thus,  we  first  focused  on  these  tasks. 

The  aim  was  to  initialize  the  spin(s)  at  the  end  of  the  chain  in  a  state  of  interest  for  the  transfer  of  either 
classical  or  quantum  information,  while  leaving  the  rest  of  the  spin  chain  in  the  maximally  mixed  state.  In 
the  first  case,  we  would  like  to  prepare  the  state  Sp\  ~  a\  0  11 ;  whereas  in  the  second  case  we  would  like  to 

prepare  one  of  the  logical  states,  e.g.,  Spy  =  cry(7x+crx<Ty  0  llri_2  [5].  Exploiting  the  spin  natural  dynamics 
and  a  combination  of  coherent  and  incoherent  control  it  is  possible  not  only  to  prepare,  but  even  to  detect 
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these  types  of  states  [5].  This  was  one  critical  step  toward  the  demonstration  of  QST  in  a  solid-state  NMR 
platform. 

The  key  insight  was  to  realize  that  even  in  the  absence  of  frequency  addressability,  the  dynamics  of  the  end- 
chain  spins  under  the  internal  dipolar  Hamiltonian  is  different  from  the  bulk  spins,  as  the  end-spins  have 
only  one  nearest  neighbor. 

Polarization  initially  in  the  transverse  plane,  5p  =  J2k= 1  ax  (prepared  from  the  thermal  state  by  a  n/2 
pulse),  evolves  under  the  internal  dipolar  Hamiltonian  at  different  rates.  The  end-spin  evolution  rate  is 
slower  by  a  factor  «  1/  y/2  as  compared  to  the  rest  of  the  chain,  due  to  fewer  numbers  of  couplings  with 
neighboring  spins.  Thus,  there  exist  a  time  t\  when  the  state  of  the  end- spins  is  still  mainly  <jx ,  whereas 
the  rest  of  the  spins  have  evolved  to  many-body  correlations.  A  second  7t/2  pulse  brings  the  end-spin 
magnetization  back  to  the  longitudinal  axis,  while  an  appropriate  phase  cycling  scheme  cancels  out  other 
terms,  thus  obtaining  the  state 

$Pend  =  &p\  +  $Pz-  (16) 

We  note  that  the  chain  geometry  prevents  breaking  the  symmetry  between  spin  1  and  N.  Here  the  phase  cy¬ 
cling  achieves  a  similar  result  of  temporal  averaging  in  the  preparation  of  pseudo-pure  states.  The  sequence 
that  prepares  this  state  can  thus  be  written  as 


7T 

2 


-tl- 

cx 


(SI) 


with  a  =  {— x,  y},  to  average  out  terms  that  do  not  commute  with  the  total  magnetization  Ez  =  az-  As 
the  phase  cycling  does  not  cancel  zero-quantum  coherences,  they  will  be  the  main  source  of  errors  in  the 
initialization  scheme. 

A  similar  control  strategy  can  be  as  well  used  to  read  out  the  spins  at  the  end  of  the  chain  even  if  the 
observable  in  inductively  measured  NMR  is  the  collective  magnetization  of  the  spin  ensemble,  Ez. 
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Fig.  5:  Transport  under  the  DQ  Hamiltonian  [5].  Data  points  are  the  experimental  data,  with  error  bars  obtained 
from  the  offset  of  the  signal  from  zero.The  lines  are  the  fits  using  the  analytical  model.  The  fitting.  LEFT:  Blue: 
Initial  state  5pth\  readout,  collective  magnetization,  T,z.  Red:  Initial  state,  5pend\  readout,  end  readout.  The  two 
curves  highlight  the  differences  arising  from  the  different  initial  state  and  readouts.  RIGHT:  Blue,  open  circles:  Initial 
state,  5pend ;  readout,  collective  magnetization,  T,z.  Red,  filled  circles:  Initial  state,  8pth\  readout,  end-spin  readout. 
The  experimental  data  shows  remarkable  agreement  between  the  two  schemes,  thus  confirming  the  validity  of  the 
initialization  and  readout  methods. 


To  measure  a  different  observable,  the  desired  state  must  be  prepared  prior  to  acquisition.  Thus  we  want 
to  turn  into  the  end-chain  state,  Eq.  (16).  In  general,  the  sequence  used  for  readout  cannot  be  a  simple 
inversion  of  the  end-selection  step  since  this  is  not  a  unitary  -reversible-  operation.  It  is  however  sufficient 
to  ensure  that  the  state  prior  to  the  end- selection  sequence  has  contributions  mainly  from  population  terms 
(oc  cr*?)  for  the  sequence  (SI)  to  work  as  a  readout  step.  Since  the  states  created  by  evolution  under  the  DQ 
Hamiltonian  are  already  of  the  form  OC  (Tkz,  the  (SI)  sequence  with  a  two-step  phase  cycling  is  enough  for 
the  end-readout  step. 
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We  verified  the  effectiveness  of  the  initialization  and  readout  methods  by  probing  the  transport  dynamics 
driven  by  the  transport  Hamiltonian,  comparing  the  end-polarized  states  and  observables  with  the  thermal- 
equilibrium  state.  In  the  timescale  of  the  experiment,  the  initial  perturbation  travels  across  ^17  spins, 
however  only  polarization  leaving  one  end  of  the  chain  could  be  observed:  a  clear  signature  of  the  polariza¬ 
tion  reaching  the  other  end  is  erased  by  the  distribution  of  chain  lengths.  Still,  the  experimental  verification 
of  initial  state  preparation  is  possible  even  at  these  short  time  scales  thanks  to  marked  differences  in  the 
signal  arising  from  the  evolution  of  thermal  and  end-polarized  states  under  DQ  Hamiltonian. 

The  left  panel  of  figure  (5)  (blue)  shows  the  observed  evolution  of  the  collective  magnetization  E*  under 
the  DQ  Hamiltonian,  starting  from  the  thermal  initial  state,  Spth  =  Modeling  the  physical  spin  system 
by  an  ensemble  of  equivalent  and  independent  spin  chains  with  nearest-neighbour  couplings  only,  we  can 
derive  analytical  formulas  for  the  evolution  to  fit  the  experimental  data,  finding  a  very  good  agreement.  The 
red  data  (5)  show  instead  the  evolution  of  the  end  polarized  initial  state  under  the  transport  Hamiltonian, 
measured  using  the  readout  strategies  outlined  above.  This  experiment  is  thus  a  direct  simulation  of  quantum 
state  transport.  The  two  sets  of  data  show  very  different  chain  dynamics  for  the  two  initial  states  (with  and 
without  end  selection),  giving  an  experimental  validation  of  our  initialization  method. 

To  further  validate  the  initialization  and  readout  method,  in  the  right  panel  of  figure  (5)  (blue,  open  circles) 
we  plot  the  system  dynamics  when  starting  from  an  end-polarized  state  (where  polarization  is  localized  at 
the  ends  of  the  chain)  and  reading  out  the  collective  magnetization.  The  red  (filled  circle)  data  shows  a 
complementary  measurement  where  we  start  from  thermal  initial  state,  given  by  the  collective  magnetiza¬ 
tion,  and  read  out  the  ends  of  the  chains  after  evolution  under  the  DQ  Hamiltonian.  Both  these  data  sets 
were  fitted  by  the  analytical  expression  that  can  be  found  in  the  nearest-neighbor  coupling  limit,  thanks 
to  a  Jordan- Wigner  transformation  to  non-interacting  fermions.  Furthermore,  the  data  and  fittings  for  end 
selection  and  end  readout  measurements  are  very  similar.  This  indicates  the  robustness  of  the  readout  step. 
We  use  the  end-selection  scheme  not  only  to  prepare  the  end-chain  polarized  states,  but  also  the  logical  states 
required  for  quantum  information  transport.  First  we  prepare  the  end  polarized  state  5  pend  by  the  sequence 
(SI).  Then  this  evolves  under  the  DQ  Hamiltonian  for  a  very  short  time  £dq  =  14.7/iS,  thus  creating  a 
two-spin  correlated  state  as  required.  We  can  write  this  initialization  sequence  as 
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—  h  — 

a 
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—  dq7, 
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(S2) 


where  DQ7_X  is  propagation  under  Umq ( t dq  )  =  e  lU<>Qfi'>Q  and  DQ7-'7  under  u\4q.  Setting  [a,/3, 7]  = 
[— x ,  x ,  x\,  the  state  after  the  sequence  (S2)  is  approximately  given  by  zero  and  double  quantum  coherences, 
Spend(.tm)  ~  cr^2  +  an-l,n  +  ^2  +  an-l,n’  where  =  4x4  “  4°*)  and  °ij  =  0.x 4  +  0,4). 

A  double  quantum  filter  given  by  the  four-step  phase  cycling  scheme, 


[a,/3,7]  =  {[-x,x,x\]  [y,-y,x\;  [~x,~x,y]-,  [y,y,y]} 

cancels  out  the  zero-quantum  terms  and  selects  the  double-quantum  terms,  which  is  our  desired  state:  Spy  oc 
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Fig.  6:  Evolution  of  the  logical  state  5py  under  the  transport  Hamiltonian,  first  reported  in  [5].  The  logical  initial 
state  was  prepared  using  the  sequence  (S2)  and  its  evolution  under  the  DQ  Hamiltonian  monitored  by  observing  the 
collective  magnetization.  Circles  represent  the  experimental  data  and  the  dashed  line  is  the  fit  to  the  analytical  model. 
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o\  2  +  ^Q1  n.  Figure  (6)  shows  the  evolution  of  this  state  under  the  DQ  Hamiltonian  while  the  dynamics 
was  monitored  by  measuring  the  collective  magnetization,  SL  oc  Tr  ^UMQ^PyU^MQaz^.  We  note  that  this 
experiment  implements  the  transport  of  quantum  information  via  a  maximally  mixed  quantum  channel. 

3.2  Decay  of  spin  coherences  in  one-dimensional  spin  systems 

While  working  to  build  the  control  tools  to  realize  large  scale  quantum  information  architectures,  we  also 
investigated  their  coherence  properties  [6],  as  multi-qubit  systems  in  correlated  and  entangled  states  are 
usually  more  fragile  to  decoherence.  We  experimentally  studied  the  decay  of  such  multi-qubit  states  under 
the  action  of  a  correlated  spin  bath,  investigating  the  decay  rate  dependence  on  the  correlations  in  a  multi¬ 
qubit  spin  state.  We  leveraged  the  low  dimensionality  of  the  system  studied  -  the  linear  coupling  geometry 
provided  by  nuclear  spins  in  apatite  crystals  -  to  gain  insight  into  both  the  many-body  states  created  by  the 
coherent  Hamiltonian  dynamics  and  their  subsequent  decay.  We  were  thus  able  to  derive  a  simple  analytical 
model  that  captures  the  essential  features  of  the  multi-qubit  decays,  and  compares  well  with  the  experimental 
data.  These  results  will  be  helpful  in  paving  the  way  for  the  future  design  of  schemes  to  mitigate  the  decay. 
The  experimental  scheme  is  shown  in  Fig.  7.  The  system  is  first  prepared  in  a  suitable  initial  state  pi. 
Evolution  under  a  propagator  ?7mq  for  a  time  r  creates  a  complex,  multiple-quantum  coherence  state.  The 
system  is  then  let  evolve  freely  for  a  time  t ,  during  which  the  coherences  decay  mainly  under  the  effects  of 
the  dipolar  Hamiltonian.  In  order  to  observe  this  decay,  we  first  refocus  the  remaining  coherences  with  a 
propagator  C/^Q  before  measuring  the  spin  magnetization  via  the  usual  free  induction  decay. 

We  created  spin  correlations  by  evolution  under  the  double  quantum  (DQ)  Hamiltonian 

UdQ  =  “  o'yOy )»  (1?) 

ij 

which  is  known  to  generate  quantum  coherences  among  the  spins.  The  density  operator  created  by  evolution 
under  the  DQ  Hamiltonian  can  be  decomposed  into  its  multiple  quantum  coherences  (MQC)  components, 

p(r)  =  Umq^p^U^t)  =  (18) 

m 

where  a  multiple  quantum  term  of  order  m,  p ,  acquires  a  phase  rrup  under  a  collective  Ez  rotation  by  an 
angle  p.  The  correlated  spin  states  created  under  T^dq  evolution  contain  in  general  all  even  M  coherence 
orders.  However,  since  standard  NMR  techniques  measure  only  single-quantum  coherences  (SQC),  in  order 
to  probe  the  higher  spin  coherences  it  is  necessary  to  indirectly  encode  their  signatures  into  SQCs  which 
can  be  measured  inductively.  This  is  achieved  by  labeling  each  coherence  order  with  a  different  phase  v  by 
means  of  collective  rotations  =  exp(— ipTjz/2)  about  the  z  axis,  effectively  creating  the  phase  shifted 
DQ  Hamiltonian,  =  U^T-LbqU^.  Finally,  MQC  are  refocused  back  to  single-spin  single-quantum 
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Fig.  7:  Experimental  scheme.  The  system  is  first  prepared  in  an  initial  state  of  interest,  for  example  the  thermal 
equilibrium  state  or  ax  using  a  7r/2-pulse  (red  bar).  Evolution  via  the  DQ  Hamiltonian  Hbq  (obtained  by  a  multi¬ 
pulse  sequence,  blue  rectangles)  creates  spin  correlations  during  the  time  r.  A  phase  shift  p  of  the  propagator  encodes 
information  about  the  multiple  quantum  coherence  intensities  created.  The  state  undergoes  decay  under  the  dipolar 
interaction  Hdip  during  the  time  t.  The  correlated  state  is  refocused  by  the  inverse  propagator  U^q  before  a  7r/2  pulse 
is  used  to  detect  the  spin  free-induction  decay. 
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Fig.  8:  Normalized  decay  of  quantum  coherences  created  from  the  thermal  initial  state  by  evolving  for  a  time  r  under 
the  Hdq  Hamiltonian.  (ZQ  (OQ)  blue  open  circles,  DQ  (2Q)  red  dots,  Total  signal  black  stars).  Data  points  are  the 
normalized  signal  intensities  for  evolution  time  48//S  (left)  and  589/iS  (right).  The  dashed  lines  are  fitting  using  the 
function  in  Eq.  20.  For  comparison  we  plot  the  free  induction  decay  (green  squares)  that  we  fit  with  the  function 


(1  —  C)sinc(m2t)  e  rn^t‘2/2  -f  C  ,  with  the  second  moment  given  by  M  =  mi  +  m\j3. 


terms  and  free  induction  decay  is  measured.  Each  measurement  is  repeated  while  incrementing  p  from  0  to 
27 r  in  steps  of  5p  =  27t/2 K  where  K  is  the  highest  order  of  MQC  we  wish  to  encode. 

Since  often  the  observable  SpQ  is  proportional  to  the  initial  state  (as  it  is  the  case  for  the  thermal  equilibrium 
state  and  the  total  magnetization  along  the  z-axis)  we  can  write  the  measured  signal  as  a  correlation  S(t,r)  = 
Tr  {5p(t,  r)5p0(r)},  between  the  state  prepared  by  the  DQ  evolution,  5p0(r)  =  Umq(t)5p0uI/[q(t)  and 

the  same  state  after  decay  under  the  dipolar  evolution,  5p(t ,  r)  =  UdipSp0{T)u\ip.  The  signal  intensities  of 
various  coherence  orders  are  given  by  the  Fourier  Transform  with  respect  to  the  phase  p: 

K 

/(m)(f,r)  =  Tr(  5pp\T)5p(m\t,T)}  =  ^2sk(t,r)e-ikm^,  (19) 

k= 1 


where  Sk(t,  r)  =  Tr  r)5p0j  is  the  signal  acquired  in  the  kth  measurement  when  setting  p  =  nk/K. 

We  studied  the  decay  of  MQC  intensities  created  under  T~Lbq  starting  from  an  initial  thermal  state  ( 5pth  ~ 
E*)  and  other  states  of  interest  (end-polarized  states  and  transverse  polarization  states).  The  decoherence 
dynamics  was  studied  by  repeating  the  experimental  scheme  described  above  while  varying  the  DQ  evolu¬ 
tion  time  (t)  from  36ps  to  925/xs  and  the  decay  time  (t)  from  0  to  145/is  (which  is  on  the  order  of  the  free 
induction  decay  time). 

We  fitted  the  decay  curves  to  Gaussian  functions, 

G(t,  r)  =  A(t)  ([1  -  C(r)}e~M^t2/2  +  C(r))  ,  (20) 

where  A  (amplitude),  M  (second  moment)  and  C  (asymptote)  are  used  as  fitting  parameters  that  vary  with 
the  DQ  time  r.  As  shown  by  the  behavior  of  the  fitting  parameters  in  Fig.  9,  the  system  exhibits  an  interesting 
dynamics  as  a  function  of  the  DQ-time  evolution.  This  is  in  contrast  to  what  was  observed  in  3D  systems, 
where  the  decay  simply  becomes  monotonically  faster  as  the  DQ-time  r  is  increased.  The  difference  can  be 
traced  to  the  fact  that  the  constrained  coupling  topology  in  ID  systems  allow  for  a  slower  decay  dominated 
by  nearest-neighbor  interactions,  while  in  3D  systems  the  decay  is  more  rapid  and  diffusion-like. 

We  were  able  to  derive  an  analytical  model  for  the  parameters  of  this  decay,  the  asymptote  C 

C(T)  =  =  ^-2y^cos(4rcosK)^  (21) 
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Fig.  9:  Left:  Asymptotes  of  the  experimental  decay  curves  (Czq  blue,  Cdq  red  Ctotai  black).  The  curves  were  fitted 
by  Eq.  21  with  an  additional  linear  term  in  r,  giving  a  dipolar  coupling  of  7.676  x  103rad/s.  Right:  Decay  rates  of  the 
different  MQ  components  as  a  function  of  double  quantum  evolution  time  r.  Points  are  experimental  results  obtained 
from  the  fitting  of  individual  decay  curves  for  each  DQ-time  r.  Error  bars  are  estimated  from  fitting  of  decay  curves 
with  20).  The  dashed  lines  are  /i  M(b,  r),  where  M(6,  r)  are  the  analytical  curves. 


and  the  decay  moment  M  =  Mzz  +  2MXX: 


Mzz  = 


16 

~N 


N 


52  I/p.<z(2t)|2  “  52  (l/<z,<z-i(2r)|2  +  l/?,?-2(2r)|2  +  |/i;?(2r)|2) 

p/g  q= 2 


(22) 


MXx  —  jy  jy  ^  [fp+i,q(2r)  /p-i,g(2r)]  [/p,g+i(2r)  /p,g-i(2r)]  ,  (23) 

p,q 

with 

/OT(r)  =  ^-rE(-l)P  sin(^)sin(^)e-2iWcosK,  (24) 

k 

where  N  is  the  chain  length  and  k  =  -^L-.  These  expressions  were  obtained  by  approximating  the  Hamil¬ 
tonian  with  a  nearest-neighbor  only  interaction  and  using  a  mapping  to  non-interacting  fermions. 

We  further  investigated  the  dependence  of  decoherence  rate  for  different  initial  states  that  evolve  under  the 
DQ-Hamiltonian  to  interesting  many-spin  states  showing  different  types  of  spin-spin  correlations.  Since  the 
decay  process  is  non-Markovian,  but  it  is  due  to  a  highly  correlated  spin  bath,  we  found  a  very  rich  dynamics, 
where  decoherence  rates  (quantified  by  the  second  moment  of  the  the  decay)  depend  in  a  non-trivial  way  on 
the  degree  of  localization  of  the  state  as  well  as  on  its  coherence  with  respect  to  the  quantization  basis.  In 
particular  we  found  that  large  spin  clusters,  with  correlations  established  among  many  spins,  decay  faster 
under  a  correlated  bath,  even  if  their  coherence  order  is  not  very  large.  This  is  in  contrast  to  the  decay 
under  simple  dephasing,  where  the  coherence  order  (and  for  pure  states,  the  entanglement)  is  critical  in 
determining  the  decay  rate.  While  it  was  not  possible  to  separate  the  coherence  order  and  the  number  of 
correlated  spins  in  the  dynamics  of  3D  spin  systems  (as  they  grow  at  the  same  time),  here  we  were  able 
to  get  more  insight  by  using  spin  chains  and  exploring  different  initial  states.  In  addition  we  found  that 
restricting  the  dynamics  in  one  dimension  slows  down  the  decay,  which  could  be  beneficial  to  create  larger 
coherent  quantum  states. 


3.3  Engineering  Dephasing  and  Depolarizing  Noise 

In  the  first  year  of  the  project,  we  focused  most  of  our  efforts  in  devising  strategies  to  achieve  perfect 
quantum  information  transport  by  engineering  the  spin  interactions.  We  found  that  dynamically  modulating 
the  spin-spin  couplings  could  lead  to  perfect  fidelity  of  information  transfer  even  in  disordered  systems. 
However,  these  results  could  be  voided  if  decoherence  is  at  play  [33].  Thus,  we  later  focused  our  attention 
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Fig.  10:  Unit  cell  of  the  fluorapatite  crystal  [CasCPO^F],  high¬ 
lighting  the  geometry  of  the  fluorine  chains  (in  red)  and  the  phos¬ 
phorus  nuclei  (in  yellow).  Each  fluorine  spin  is  surrounded  by  3 
phosphorus  spins.  As  the  phosphorus  are  in  a  thermal  state,  they 
create  a  local  random  magnetic  field  at  the  location  of  each  flu¬ 
orine  spin  that  induces  decoherence  and  localization  unless  the 
spins  are  decoupled  by  active  control. 


to  studying  effects  of  decoherence  and  disorder  during  the  transport.  To  do  so,  we  exploited  the  presence  of 
a  second  nuclear  species  in  our  experimental  system. 


Fluorapatite  crystals  indeed  contain  phosphorus,  whose  nucleus  31 P  has  a  spin- 1/2  with  gyromagnetic  ratio 
7 p  —  17.25MHz/T  (compare  to  the  fluorine’s  7 p  —  40MHz/T).  While  in  previous  experiments  the  system 
dynamics  was  protected  against  the  fluorapatite  31 P  spin  bath,  it  is  interesting  to  study  what  would  be  its 
effects  on  the  transport  dynamics.  In  particular,  disorder  and  noise  have  been  indicated  as  potential  sources 
of  localization,  which  would  prevent  the  goal  of  transferring  quantum  information.  The  couplings  to  31 P  are 
weaker  than  F-F  in-chain  couplings,  since  the  gyromagnetic  ratio  is  smaller  and  distance  is  larger  (the  closest 
distance  between  a  19F  and  a  31P  is  3.5A  (compare  to  a  distance  3.4A  between  Fluorines  in  the  spin  chains). 
However,  since  there  are  3  31 P  per  each  19F,  their  combined  effect  is  large,  as  indicated  by  calculations  of 
the  second  moments.  The  31 P  nuclei  act  as  a  spin  bath,  coupling  to  the  fluorine  via  the  Hamiltonian 


upf  =  53 

ji 


3cOs(l?j;)2  -  1  p  F 
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From  the  point  of  view  of  the  fluorine  spins,  this  amount  to  the  presence  of  a  local  magnetic  field 

Bj  7P7F%0  3  008(0 ji)2  -  1  !  p  \ 

p  4^  Vji\^  ' 

which  depends  on  the  (thermal)  state  of  the  phosphorus  spins.  In  addition  the  31 P  interact  among  themselves, 
thus  giving  rise  to  a  time-dependent  magnetic  field.  Using  multiple  pulse  sequences,  we  can  engineer  the 
19F  dynamics  (for  example  engineering  the  transport  Hamiltonian  Hdq  while  at  the  same  time  refocusing 
or  shaping  the  spin  bath  interaction. 

This  allows  us  to  study  not  only  decoherence,  but  more  importantly,  the  effects  of  a  disorder  onto  the  transfer 
of  information.  In  particular,  we  are  able  to  observe  experimentally  Anderson  localization. 


Fig.  11:  Simulation  of  localization  in 
spin  chains  due  to  a  disordered  mag¬ 
netic  field.  Here  we  plot  a  measur¬ 
able  observable,  the  many-body  cor¬ 
relation  length,  describing  the  average 
number  of  correlated  spins.  The  cor¬ 
relation  length  can  be  measured  ex¬ 
perimentally  from  the  multiple  quan¬ 
tum  coherence  intensities. 


3.4  Polarization  Generation  and  Transfer  in  Electronic  spin  systems 

In  addition  to  nuclear  spin  chains  in  fluorapatite,  we  investigate  an  alternative  spin  system,  connected  to  the 
NV  center  in  diamond.  While  in  diamond  it  is  not  yet  possible  to  obtain  regular  arrays  of  spins,  we  still 
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Fig.  12:  Left:  Hartmann-Hahn  matching  sequence  for  polarization  contact.  Right:  Energy  matching  in  the  rotating 
(dressed)  frame,  here  for  two  spin- 1/2.  Measured  decay  of  NV  spin-lock  signal  as  a  function  of  spin-lock  duration. 
When  the  PI  bath  is  not  driven  we  record  the  blue  trace.  Driving  the  PI  bath  simultaneously  at  all  5  ESR  resonances 
such  that  the  collective  PI  bath  Rabi  frequency  equals  the  NV  Rabi  frequency  (8  MHz)  gives  the  red  trace,  indicating 
strong  NV/P1  polarization  transfer.  Solid  lines  represent  fits  to  decaying  exponentials. 


started  investigating  polarization  transfer.  This  work  laid  the  foundations  for  further  experiments  that  will 
combine  our  results  in  control  of  local  quantum  registers  and  results  in  quantum  information  transfer  with 
spin  wires. 

We  investigated  generation  of  polarization  and  its  transport  by  exploiting  the  NV  center  in  diamond  as  a 
source  of  polarization  and  the  bath  of  electronic  spins  as  the  bus.  The  most  common  nitrogen  defect  in 
diamond  is  the  PI  center,  a  nitrogen  substitutional  atom  with  an  unpaired  electron.  The  PI  center  has 
electronic  spin- 1/2,  with  a  strong  hyperfine  interaction  with  its  nuclear  spin  giving  rise  to  5  distinct  resonance 
frequencies.  Polarization  transfer  between  the  NV  center  and  other  spins  can  be  achieved  by  continuous 
irradiation  at  the  Hartmann-Hahn  matching  condition  [34].  The  CW  field  effectively  establishes  a  resonance 
condition  in  the  rotating  frame  (or  dressed  basis). 

As  the  NV  spin  can  be  polarized  optically  with  high  efficiency,  the  polarization  can  be  transferred  to  the 
electronic  spin  bath.  The  transfer  was  implemented  by  achieving  an  energy  matching  between  the  NV 
and  PI  spins  in  the  rotating  frame.  We  proved  this  mechanism  experimentally,  overcoming  the  difficulties 
associated  with  the  large  hyperfine  of  the  PI  centers,  which  required  us  to  simultaneously  drive  5  resonance 
transitions  [35].  The  energy  contact  was  confirmed  by  a  reduction  in  the  NV  coherence  time  (Fig.  12). 
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Abstract 

The  goal  of  the  project  "Modular  Paradigm  for  Scalable  Quantum  Information"  was  to  address  some  of  the 
challenges  facing  the  field  of  quantum  information  science  (QIS).  The  requirements  to  reliably  control  a 
scalable  quantum  system  while  staving  off  decoherence  pose  a  contradiction,  as  fast  control  implies  a 
strong  coupling  to  a  controlling  (external)  system,  but  this  entails  an  undesired  interaction  with  the 
environment,  leading  to  decoherence.  We  thus  studied  a  new  paradigm  for  QIS,  given  by  smaller  building 
blocks  (composed  of  a  quantum  register  and  an  actuator)  that  are  connected  by  spin  wires. 

The  project  had  activities  along  these  two  lines  of  research. 

On  one  side,  we  studied  the  working  principles  of  the  actuator/register  model,  in  particular  its  control  and 
the  indirect  control  of  the  quantum  register. 

We  developed  a  comprehensive  theory  of  time-optimal  control  of  the  quantum  register,  finding  a  minimal 
set  of  parameters  that  need  to  be  optimized  to  obtain  the  time-optimal  control  sequence  for  any  desired 
quantum  gate.  This  enables  us  to  devise  an  efficient  numerical  method  to  find  the  time-optimal  sequence 
and  in  turn  to  investigate  its  performance  on  a  realistic  system. 

In  addition,  we  developed  experimental  control  of  an  exemplary  actuator/register  system,  the  nitrogen- 
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vacancy  (NV)  electronic  spin  center  in  diamond  and  its  associated  nitrogen  nuclear  spin. 

We  used  the  control  to  characterize  the  Hamiltonian  of  the  system  and  to  achieve  enhanced  control  speed, 
exploiting  the  presence  of  the  quantum  actuator.  In  addition,  we  showed  how  to  use  a  quantum  controller  to 
perform  quantum  feedback,  to  protect  the  qubit  systems  from  decoherence. 

The  main  results  of  this  first  task  are: 

-  We  found  the  shape  of  the  time-optimal  control  strategy  to  engineer  desired  qubit  gates  via  the  actuator 
control.  Using  algebraic  methods  we  were  able  to  find  the  time-optimal  control  solutions  for  the  most 
general  SU(2)  synthesis  problem. 

-  We  developed  an  efficient  numerical  method  to  find  the  time-optimal  control  sequence  and  applied  it  to 
find  gates  for  the  coupled  NV-nuclear  spin  system. 

-  We  exploited  an  hybrid  control  strategy,  combining  classical  and  quantum  controllers,  to  achieve  faster 
control  of  a  nuclear  spin  qubit. 

-We  exploited  a  long-lived  spin  qubit  as  a  controller  to  implement  quantum  feedback  to  preserve  the 
coherence  of  the  NV  electronic  spin 

Our  activities  demonstrated  that  quantum  controllers  can  provide  powerful  advantages  over  classical 
controllers,  achieving  faster  and  more  robust  control  and  make  new  tasks,  such  as  feedback  control, 
possible. 

For  the  second  activity  of  the  program  we  explored  quantum  information  transfer  using  quantum  spin  wires, 
which  transport  information  via  the  free  evolution  under  their  mutual  interaction. 

We  combined  a  theoretical  study  of  the  conditions  for  perfect  quantum  state  transfer  with  an  experimental 
investigation  the  role  of  decoherence  and  disorder  in  the  transport  process,  using  an  experimental 
implementation  with  nuclear  magnetic  resonance  techniques  applied  to  apatite  crystals. 

Results  we  achieved  pertinent  to  this  second  task  include: 

-  We  performed  the  first  experimental  implementation  of  quantum  information  transfer  in  mixed-state  spin 
chains. 

-  We  devised  protocols  for  Hamiltonian  engineering  that  can  lead  to  perfect  quantum  information  transfer 
even  in  the  absence  of  single-spin  addressability.  These  techniques  can  be  more  generally  used  for 
quantum  simulations  in  a  variety  of  quantum  systems. 

-  We  found  that  the  restriction  of  the  spin  system  to  one-dimensional  geometries  yields  a  complex,  non- 
Markovian  decoherent  dynamics.  However,  it  also  makes  the  coherence  times  longer  that  in  3D  systems, 
thus  pointing  to  advantages  to  be  found  in  particular  geometries  for  larger  quantum  information 
architectures. 


-  We  studied  the  polarization  transfer  from  NV  centers  in  diamond  (that  can  be  polarized  optically)  to  dark 
spins  associated  with  single  Nitrogen  defects  (PI  centers)  in  diamond.  We  positively  assessed  the 
potential  of  this  system  for  a  scalable  quantum  information  architecture. 

Our  results  provide  new  powerful  tools  for  achieving  robust  quantum  wires  and  strategies  for  counteracting 
the  effects  of  noise  and  disorder.  At  the  same  time,  we  devised  novel  techniques  that  have  broader  impacts 
on  quantum  simulations,  from  Hamiltonian  engineering  to  novel  ways  of  engineering  the  bath. 
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